Ferromagnetic behavior in magnetized plasmas 
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We consider a low-temperature plasma within a newly developed MHD Fluid model. In addition 
to the standard terms, the electron spin, quantum particle dispersion and degeneracy effects are 
included. It turns out that the electron spin properties can give rise to Ferromagnetic behavior in 
certain regimes. If additional conditions are fulfilled, a homogenous magnetized plasma can even be 
unstable. This happen in the low-temperature high-density regime, when the magnetic properties 
associated with the spin can overcome the stabilizing effects of the thermal and Fermi pressure, to 
cause a Jeans like instability. 

PACS numbers: 52.27.-h, 52.27.Gr, 67.57.Lm 



Already in the 1960's, Pines studied the excitation 
spectrum of quantum plasmas for which we have a 
high density and a low temperature as compared to nor- 
mal plasmas. In such systems, the finite width of the 
electron wave function makes quantum tunnelling effects 
crucial, leading to an altered dispersion relation. Since 
the pioneering paper by Pines, a number or theoretical 
studies of quantum statistical properties of plasmas has 
been done (see e.g. Ref. and references therein). 
For example, Bezzcrides & DuBois presented a kinetic 
theory for the quantum electrodynamical properties of 
nonthermal plasmas 0], while Hakim & Heyvaerts pre- 
sented a covariant Wigner function approach for rela- 
tivistic quantum plasmas Recently there has been 
an increased interest in the properties of quantum plas- 
mas, e.g. 0, S, 0, S, H, 0, [III, El- The studies has been 
motivated by the recent development in nanostructured 



the discovery of 




materials [l3 and quantum wells 

ultracold plasmas [l5| (see Ref. [l6| for an experimental 
demonstration of quantum plasma oscillations in Ryd- 
berg systems), astrophysical applications or a gen- 
eral theoretical interest. Moreover, it has recently been 
experimentally shown that quantum dispersive effects are 
important in inertial confinement plasmas [l8| . The list 
of quantum mechanical effects that can be included in 
a fluid picture includes the dispersive particle proper- 
ties accounted for by the Bohm potential [a, [gI [3, [sj, |9| , 
the zero temperature Fermi pressure |M[M13,Ml3|, spin 
as well as certain quantum electro- 
20l I2IL Within such descriptions 
nil . I19I . |20|, [21 1 , quantum and clas- 
sical collective effects can be described within a unified 
picture. 

In the present paper we will make use of general equa- 
tions for spin plasmas that were derived in Ref. llOll , and 
further developed towards the MHD regime in [ll||. Ex- 
ploring the basic set of equations presented in Ref. [llj we 
demonstrate that the standard plasma behavior can be 
significantly changed by the electron spin properties, and 
that the plasma can even show Ferromagnetic behavior 
in the low temperature high density regime. Further- 
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more, a homogeneous magnetized plasma can actually 
be unstable, also when the spin degree of freedom is in 
thermodynamic equilibrium. The instability is due to 
the magnetic attraction of spins, and the mechanism is 
conceptually similar to the well-known Jeans instability 
[2^ . Applications of our results to laboratory and astro- 
physical plasmas are discussed. 

Adopting the spin-MHD equations put forward in Ref. 
[m . our plasma is described by the continuity equation 



dp 

dt 



+ V • (fiv) = 0, 
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the momentum equation 
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and the idealized Ohm's law 
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where p is the plasma density, v the fluid velocity, B 
the magnetic field, M the magnetization, P the pres- 
sure, me, (mi) denotes the electron (ion) mass and h 
is Planck's constant. In addition to the standard ideal 
MHD momentum equation, Eq. ^ contains the quan- 
tum Bohm potential (which tends to smooth the density 
profile), as well as magnetization effects due to the elec- 
tron spin. Equations (U)-® should be complemented by 
an expression for the Magnetization, as well as an equa- 
tion of state for the pressure. In thermodynamic equi- 
librium, the degree of spin alignment with an external 
magnetic field is described by the Brillouin functions Bs, 
where the index s is the spin number. For spin-^ parti- 
cles we have Bi^^ipsB /T) = tanh(/iB_B/T), leading to 
a corresponding Magnetization 



M 



^ tanh ( B. 



mi 



V T 



(4) 



2 



Here B denotes the magnitude of the magnetic field and 
B is a unit vector in the direction of the magnetic field, 
fiB = eft/2me is the Bohr magneton, e is the magni- 
tude of the elementary charge, and T is the temperature 
given in energy units. In general the argument of the 
tanh-function can vary, if for example the magnetic field 
strength varies. However, in case the variations of the 
magnetic field occurs on a time-scale shorter than the 
characteristic spin relaxation time, the degree of align- 
ment can be considered as approximately constant. Since 
spontaneous spin changes does not occur for single elec- 
trons (due to angular momentum conservation) , this spin 
relaxation time is not smaller than the inverse collision 
frequency, which can be taken as infinite in many appli- 
cations. This case will be considered for the reminder of 
this article, and thus tanh(yUBi?/T) tanh(/iBi?o/7o) 
in Eq. ([4]), where the index denotes the unperturbed 
background value. Furthermore, for low electron temper- 
atures, it is necessary to include the contribution from 
the zero temperature Fermi pressure in the equation of 
state. Writing the equation of state as 



VP = ciVn 



(5) 



we emphasize that the ion-acoustic velocity Cg includes 
the contribution from the ion and electron thermal mo- 
tion, as well as the electron Fermi pressure. Thus we 
have 



(6) 



where vu and Vte are the (effective) ion and electron ther- 
mal velocities 23|, whereas vpe is the electron Fermi ve- 
locity [2J|. Equations H])-® constitute a closed set that 
describe the spin modified quantum MHD equations. 

In what follows, we will study the linear modes of this 
system, with a particular focus on the stability proper- 
ties. With n = no + ni, B = Bq + Bi, M = Mq + Mi, 
and V = Vi, such that ni <C uq, \Bi\ <C |-Bi|, and 
\Mi\ ^ 1-^0 1 7 we linearize our equations in the per- 
turbed variables. Assuming that the background quan- 
tities are constants, the general dispersion relation can, 
after a Fourier decomposition, be written 
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where Ca is the spin-modified Alfvcn velocity given by 

Ca 



Ca = 



1/2 



, (8) 



1 + {hw^Jmc-^Loie^)tSinh{fiBBo/To) 

Ca is the standard Alfvcn velocity Ca ~ {B'n/ fJ-oPoY^^ 

Mb -Bo 
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and 



Here 
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(nge^/eo™)^^^ is the plasma frequency. 



,(0) 



is the electron cyclotron frequency associated with the 
external magnetic field (i.e. with the contribution to Bq 
from the spin sources excluded). The relation between 
the full electron cyclotron frequency LOce = eBo/me and 
Wce^ is given by = Wc? + hiL!pptanh{jj,BBo/To)/mc'^. 
We stress that Vsa, which to some extent can be consid- 
ered as an effective acoustic velocity, may be imaginary 
for a strongly magnetized plasma due to the spin con- 
tribution, a fact which will be explored in some detail 
below. 

In deducing Eq. ([7]) we have assumed that the spin 
orientation has reached the thermodynamic equilibrium 
state in response to the external magnetic field. This en- 
sures there is no free energy stored into the spin degree 
of freedom, and as a consequence it turns out that the 
shear Alfven mode described by the first factor of ([7]) is 
always stable, since clearly Ca is always real. This is re- 
lated to the fact that this particular mode has no density 
perturbations. By contrast the second factor, describing 
the fast and slow magnetosonic modes, does not neces- 
sarily predict stability. The reason is that the electrons 
carry spin, and that they thus behave as single magnets 
to some extent. Just like magnets or gravitating mat- 
ter, the electrons may thus attract each other, leading to 
an exponentially growing density, similar to the gravita- 
tional Jeans instability. Naturally electrostatic repulsion 
among the electrons could in principle act as a strong 
counteracting force to this scenario. However, within the 
low- frequency MHD limit, ions and electrons move to- 
gether, and thus the Coulomb force does not provide a 
stabilizing mechanism. To shed some further light on the 
stability properties, we consider propagation perpendic- 
ular to the external magnetic field, which is the geome- 
try which leads to instability most easily. For the case 
k = kx, Eq. ^ reduces to 
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The condition for instability is thus that the last nega- 
tive term of pT|) dominates over all the others. Under 
this assumption, we have depicted the growth rate as a 
function of k in Fig. 1. 

The necessary and sufficient instability condition can 
thus be written as 



Psp + Pm + P+ ^^^^ < 0, 

4me 



(12) 
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FIG. 1: The growth rate Imtj as a function of k obtained 
from the dispersion relation (11). The relation (11) is of 
the form u = k{h^k^ — a^)^'^^, and we use the normalization 
uj = |b|Lj/|a|^ and A; = |&|A:/|a|. We note that we have a max- 
imal growth rate for the wavenumber fcmax = |a|/(v^|b|) ~ 
[2me[(|Psp| — PiTL — Py\/noh']^^'^ , corresponding to a maximal 
growth rate 7max = (me/mi)^/^[|Psp| - Pm - P)]/noh. Fur- 
thermore, although (11) is a linear dispersion relation, the 
present instability shows similarities to the modulational in- 
stability. 



where the total pressure Ptot = ^sp + ^rn + 
P consists of the effective spin pressure Pgp = 
—nQhuJcei8.iih{iJ,B Bo /To) which is the only negative pres- 
sure term and therefore the source of the instabihty, 
the magnetic pressure Pm and the ordinary pressure 
P = norriiCg, containing both the thermal and Fermi 



pressure part, 
is given by 

P„, = 



Furthermore, the magnetic pressure P„ 
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Thus a necessary (although not sufficient) condition for 
instability is 



A > Ar, = 



-,1/2 
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(14) 



which means that the instability is stabilized for short 
wavelengths A = 27r/fc, similar to the Jeans instabil- 
ity. The stabilizing influence for short wavelengths stems 
from the Bohm potential. Furthermore, the partial in- 
stability condition 



^'spl >P, 



(15) 



means that a finite pressure also may lead to stabiliza- 
tion. We note from Eq. ([6]) that a low temperature 
is not necessary to fulfill this condition, since the zero- 
temperature Fermi velocity contributes to and thereby 
to P. However, for a sufficiently strong magnetic field, 



clearly ((T5|) can be fulfilled. Finally, the last subpart of 
the instability condition reads 



IPsdI > Prr 



(16) 



which means that the magnetic pressure also acts as a 
stabilizer. For a given magnetic field, this condition may 
be fulfilled for a sufficiently high density. However, in- 
creasing the density means that the Fermi velocity is in- 
creased, which may lead to a violation of (fTS)) . To simul- 
taneously fulfill (Uni) and (Uni), and thereby to fulfill p2|) . 
it is required that the second term in the denominator of 
the right side of p3|) is larger than unity, i.e. 
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Since the spin cannot contribute much to the unper- 
turbed field unless the temperature is small enough to 
allow a significant alignment, this condition in turn re- 
quires 



1 < 



nieC^To 



(18) 



For temperatures small enough to fulfill (jlSp . the spin 
contribution to the unperturbed field dominates over the 
external field, and the plasma thus shows Ferromagnetic 
behavior. Contrary to a normal ferromagnet, however, 
the density variations are not restricted, which render 
possible the instability discussed above. However, plas- 
mas with the required background parameters are not 
easy to produce, as we can see from the following ex- 
amples: Firstly, if we chose a high density plasma like in 
incrtial fusion experiments, p ^ lO^kg/m'^, ferromagnetic 
behavior occurs for temperatures T < 10^ — lO^K, as de- 
scribed by the inequality p8|) . In this regime the Alfven 
velocity can differ much from the standard Alfven veloc- 
ity, as given by ([5]) . If, in addition, we want the Jeans- like 
instability to occur, the most severe condition to fulfill is 
which require temperatures T < 20K for standard 
laboratory field strengths. Until a few years ago, the only 
known plasmas where such low temperatures could be 
found were solid state plasmas, which do not fit into the 
MHD-likc model used here. However, recently gaseous 
plasmas with ultra-low temperatures, T < 10~^K, has 
been constructed with the aid of Rydberg atoms [13, [IB] ■ 
Unfortunately, the combined requirement of a reasonably 
high-density, Eq. ^5]) . rules out the spin-instabihty de- 
scribed above in such a laboratory setting. 

In addition to laboratory applications, the theories de- 
scribed above could be adopted for astrophysical pur- 
poses [l3|. In magnetar atmospheres, the strong mag- 
netic field makes it possible to fulfill the conditions 
even for a relativistic temperature. In that case we should 
adopt the theory to a pair plasma I4I rather than an ion- 
electron plasma. Furthermore, for white dwarf stars, the 
high density makes the condition fulfilled. 
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